Shielding of a moving test charge in a quantum plasma 
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The linearized potential of a moving test charge in a one-component fully degenerate fermion 
plasma is studied using the Lindhard dielectric function. The motion is found to greatly enhance 
the Friedel oscillations behind the charge, especially for velocities larger than a half of the Fermi 
velocity, in which case the asymptotic behavior of their amplitude changes from 1/r 3 to 1/r 2 ' 5 . In 
the absence of the quantum recoil (tunneling) the potential reduces to a form similar to that in a 
classical Maxwellian plasma, with a difference being that the plasma oscillations behind the charge 
at velocities larger than the Fermi velocity are not Landau-damped. 

PACS numbers: 52.25.Mq, 52.35.Fp, 03.75.Ss 



I. INTRODUCTION 

The Debye shielding of a moving test charge in a clas- 
sical plasma is one of the most fundamental problems in 
plasma physics (TT-fT5] . In the collisionless case the motion 
is known to result in the l/r 3 -dependence of the poten- 
tial at large distances [I] . The angular dependence of this 
asymptotic form is determined by the velocity distribu- 
tions of the plasma components [3], and, for instance, for 
a one-component Maxwellian plasma the potential is re- 
pulsive (for a particle of like charge) in front of the test 
charge and attractive behind it and perpendicular to the 
motion [5J [TUl HH] . In addition, behind the charge a large 
number of potential minima is formed at substantially 
suprathermal velocities [TJ HUJ US] ■ The presence of col- 
lisions can lead to the l/r 2 -dependence of the potential 
[5] , while the excitation of ion-sound waves can result in 
the formation of an oscillatory wake structure inside the 
corresponding Mach cone [T2HT4"] . 

While the above studies deal with classical plasmas, 
recently there has been a rapidly growing interest in 
quantum plasmas, motivated primarily by the develop- 
ment of nanostructured metallic and semiconductor ma- 
terials [T7H23]- Much attention has been given to a 
one-component weakly coupled fully degenerate fermion 
plasma [2"5T - |3"5] . Such a plasma is often described by the 
Lindhard dielectric function [30, 36 where the only quan- 
tum effects included are the degeneracy and the quantum 
recoil (tunneling). It can be derived using the Wigner- 
Poisson system [THl [Ml I3TH4T] , a quantum analog of the 
Vlasov-Poisson system. 

We present an investigation of the shielding of a mov- 
ing classical charge in a quantum plasma using the Lind- 
hard dielectric function. The free parameters in this 
model are the velocity of the test charge in units of the 
Fermi velocity and the plasma coupling parameter. The 
latter governs the role of the quantum recoil and should 
be small for the model to apply, as discussed in Sec. |IV[ In 
the limit of zero coupling parameter (i.e. in the absence of 
the quantum recoil) the potential is semiclassical in the 
sense that it can be found using the classical approach 
but with a degenerate velocity distribution [34]. Such a 
potential has the aforementioned l/r 3 -asymptote at large 



distances and was investigated in Ref. [42] at small ve- 
locities. At a finite coupling parameter but zero velocity 
the potential is known to have Friedel oscillations [3"0"ll43| . 
Our study extends these results to the general case of ar- 
bitrary velocity and coupling parameter and shows how 
the quantum recoil changes the semiclassical potential 
and how the motion modifies the Friedel oscillations. 



II. MODEL 

The potential around a point test charge Q moving at a 
constant velocity v in a three-dimensional plasma is given 
in the linear approximation by the formula [H 1301 144] 



tp(r) = 



exp(ik • 



_Q l_ 

47re 27r 2 J k 2 D(k-v,k) 



dk, 



(1) 



where r denotes the position relative to the instantaneous 
position of the charge, eo is the electric constant, and 
D(ui, k) is the dielectric function of the plasma. We take 
the screening to be due to the response of a single fully 
degenerate plasma component (e.g. electrons), with all 
the other components remaining fixed as a homogeneous 
neutralizing background. To describe their response, we 
use the Lindhard dielectric function [271 \M I5UI 1521 155] : 



D(w,k) 



/(p + fik/2) - /(p - fik/2) 
u + iv — k • p/m 



I .... _!L I ±in ■■,-• ;'■--' . ---' dp, 

(2) 

where w p = ^Jne 2 /(e^m) is the plasma frequency, n is the 
particle number density (of the component that responds 
to the test charge), e and m is their charge and mass, 
respectively, v is an infinitesimal positive number (i.e. 
the limit v — ► + should be taken), /(p) is the three- 
dimensional Fermi-Dirac distribution function: 



/(p) 




if |p| <p F , 
if |p| >p F , 



(3) 



Pf = h(3Tr 2 n) 1 ' 3 ' is the Fermi momentum, and fi is the 
Planck constant over 2ir. The applicability of the model 
is discussed in Sec. HVl 



Carrying out the integral in Eq. (pi) gives 



D(u/,k) = 1 
1 1 



l\i /\.rp rp ZiCL 



[F(uj + iv + a,k) - F(lo + iv - a, k)] , (4a) 



where 



2m' 



(4b) 



vp — pf /m is the Fermi velocity, and Xtf — wf/(w p v3) 
is the Thomas-Fermi screening length; the principal 
branch of the complex logarithm should be taken. 

Unless otherwise stated, the figures in this paper have 
been generated by numerical integration of Eq. (nl) using 
the dielectric function Q. We choose the free parameters 
to be the degenerate Mach number M — v/vp and the 
parameter 



'/ 



hujp 
4Ep 



(5) 



where Ep = p F /(2m) is the Fermi energy. Note that the 
parameter rj is related to the plasma coupling parameter 



- 71 / 



/(A-keqEf) via 



7lT 



'/ 



2(3tt 2 ) 2 / 3 ' 



(6) 



In Appendix [XJ we have reformulated the problem in 
non-dimensional form in terms of the free parameters rj 
and M. 



III. RESULTS 

A. Semiclassical limit 

We first consider the semiclassical limit rj — > 0, where 
the degeneracy of the unperturbed state is the only quan- 
tum effect. The extent to which the results of the present 
subsection are relevant despite the fact that we neglect 
relativistic effects is discussed in Sec. lIVI In the semiclas- 
sical limit we have 131 



D(k-v,k) = 1 



1 



K A/pp 



G(Mk-v), 



where the hat denotes unit vectors, 



nt \ 1 x + it 



x + it + 1 
x + ie — 1 



(7) 



(8) 



and e is an infinitesimal positive number. If we set v = 
we get the exponentially screened potential of the Debye 
form: 



(p(r) 



Q 



Ane^r 



exp 



\TF 



(9) 



Now, for M < 1, the reciprocal dielectric function 
1/D(k ■ v, k) does not have singularities at real k and 
hence, as shown in Ref. [1], the asymptotic potential as 

r — > oo is 



V?(r) 



Q -h{ 1 ) x2tf 



o 



1 



where 7 is the angle between r and v, and [45] 
r 2-K r \ 



(10) 



H(rf) = - 



lim 

7T Z <5^0+ ,/q 
1 



djji 



1 



G[M(p cos 7 + yr 



/i 2 sin 



0sin7)] (fJ- + iS) 3 ' 



(11) 



As M — \ 0, we can expand H (7) in powers of M, giving 
the asymptotic formula 



1 \ Q 

ip(r) = exp 

47T6or 



r 
Xtf 



Q A 



TF 



Airen r 3 



2M cos 7 - 



1 A;P(3cos 2 7- 1) 



O 



M c 



O 



M 



(12) 



as r — ► 00, M 
was given in £ 
the factor of 2 



0. The term linear in M in Eq. (12) 



, although they appear to be missing 
This asymptotic result is qualitatively 
identical (the only difference is in the numerical coeffi- 
cients in front of the two terms) to the classical case of 
a Maxwellian plasma [5]. Equation (12 1 shows that for 



small nonzero M, an attractive (for like charges) poten- 
tial forms antiparallel and perpendicular to the motion, 
whereas the potential parallel to the motion remains re- 
pulsive (but decays as 1/r 3 instead of exponentially as 
for M = 0). These features persist at velocities up to 
and including the Fermi velocity, as illustrated in Figs, [l] 
and[2j (Note that although Fig. 1 exhibits a cone-shaped 
feature in front of the charge, this is not an oscillatory 
Mach cone of the kind that occurs behind fast-moving 
charges and is shown in Fig. k3j) 

On the other hand, for M > 1, the reciprocal dielectric 
function l/D(k-v, k) has a pole for real k [the parameter 
e sets the rule for avoiding this pole in the integral (fl])]. 
Physically, this corresponds to the fact that a test charge 
moving faster than the Fermi velocity can excite plasma 
oscillations. Furthermore, since these oscillations have 
phase velocity greater than the Fermi velocity, they are 
not Landau-damped, as can be seen mathematically in 
the fact that -D(w, k) has no imaginary part at u>/k > vp- 
This is because, with the unperturbed velocity distribu- 
tion given by Eq. ([3]) , there are no particles with velocities 
greater than vp and hence no particles satisfying the res- 
onance condition to contribute to Landau damping (this 
resonance condition is, since the quantum recoil disap- 
pears in the semiclassical limit, the same as in classical 
plasma physics) [23 |32] . The result is a strong oscilla- 
tory wake (with amplitude decreasing as 1/r) behind the 
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FIG. 1. The quantity tpr 3 [in units of Q\ TF /(4:iT6o)] in the 
semiclassical case (77 — > 0), with the test charge at the origin 
and moving at speed v = 0.5vf to the right. The white 
superimposed curve denotes the boundary between positive 
and negative potential (for Q > 0, the positive potential is on 
the right of the figure). 



FIG. 3. The ratio of the potential <p(r) to the unscreened 
Coulomb potential Q/(4neor), in the semiclassical case with 
the test charge at the origin and moving at speed v — IAvf 
to the right. Also shown (dashed lines) is the "Mach cone" , 
defined by rii/|rj_| = —\/M 2 — 1. 



B. General case 
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FIG. 2. The quantity ipr 3 [in units of Q\ TF /{4-Ke )] in the 
semiclassical case, with the test charge moving at the Fermi 
velocity, in the directions parallel to the motion (f • v = 1), 
perpendicular to the motion (f ■ v = 0), and antiparallel to 
the motion (r • v = — 1). 



test charge with an infinite number of minima, as shown 
in Fig. [3| Note that an oscillatory structure behind the 
charge is also present in the case of a Maxwellian distri- 
bution at substantially suprathcrmal velocities, but it is 
Landau-damped [lOj . As a result, the number of minima 
is finite (though very large) because the l/r 3 -asymptote 
falls off slower than the exponentially damped amplitude 
of oscillations. 



1. The case of < v < vf 

As is well known |30j . at nonzero r\ the static dielec- 
tric function £>(0,k) has a non-analyticity (the "Kohn 
anomaly" [46 ) at wavenumbers |k| = 2k p, where kp = 
PF/fo is the Fermi wavenumber. The Kohn anomaly is re- 
lated to the discontinuous Fermi surface and gives rise to 
the Friedel oscillations [U| , with the potential as r — >• 00 
given by 30 



<p{r) 



QX 



TF 



36?7 4 cos(2k F r) 



47re (2 + 3T? 2 ) 2 



1 



3 ■ ( 13 ) 



For nonzero velocities, the Kohn anomaly occurs at 
wavenumbers k such that 



|k| =2/c F (l±Mv-k), 



(14) 



as illustrated in Fig. [4] The asymptotic potential as r — > 
00 is the superposition of the contribution from the Kohn 
anomaly and from the small wavenumbers (k —¥ 0); the 
latter is identical with the total semiclassical asymptotic 
potential, since the semiclassical and general forms of the 
dielectric function coincide in the limit k > 0. 

In Appendix[E] we outline the derivation of expressions 
for the asymptotic contribution resulting from the Kohn 
anomaly. The result is, for r and v parallel: 



^) = 4^ m0,A - 



1 



ZikpArprp 



-Re 



/+(!) 



exp [2k F (l + M)ir] 



1 + 2M 



-/+(o) 



1 



(15) 
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FIG. 4. The wavenumbers k at which the Kohn anomaly in 
the dielectric function D(k ■ v, k) occurs; the velocity v is 
directed horizontally. This figure is valid for any value of n 
since the axes are normalized to the Fermi wavenumber k F . 
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FIG. 5. The amplitude of oscillations of the quantity (pr 3 
[in units of QX^ f /(Aty€o)], in the limit r — > oo, behind the 
charge (green, increasing), and in front of the charge (blue, 
decreasing), as a function of test charge velocity, for values of 
7/ (in order of increasing amplitude): 0.15, 0.3, and 0.5. (The 
amplitude for r\ — 0.15 in front of the charge is too small to 
be visible on this graph). Inset: The amplitude of oscillations 
of the quantity <pr 2,5 [in units of QX^ F /(4:-Keo)] behind the 
charge, for v/v F > 0.5. 



anomaly, with 

/±00 
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2k F (l±M\n\ 



(16a) 
(16b) 



For r and v antiparallel, we get a similar expression as 
long asM< 1/2: 



ip(r) = 



Q 



4ireor 



2kp\rpp 



Re 



, ^H(tt)X 2 tf 

cxp[-2fc F (l -M)ir] _ 
J { ) 1-2M [ ' 



1 



(17) 



The difference between the two directions arises because, 
for r and v parallel, only the outer surface of the Kohn 
anomaly, i.e. |k| = 2k F (\ + M|k ■ v|), contributes to the 
asymptotic form; whereas, for r and v antiparallel, only 
the inner surface, i.e. |k| = 2k F (l — M|k-v|) contributes. 
Finally, for r and v antiparallel and 1/2 < M < 1, a 
different asymptotic form emerges, with the amplitude 
of oscillations now decaying as 1/r 25 rather than 1/r 3 : 



ip(r) = 



1 



M2X^ F VkFr 2 - 5 



x Re 



(1-*)/- 



\2M 



.k F r\ 
•■ y ^- l 2M 



1 



(18) 



This is related to the change in concavity of the inner 
surface of the Kohn anomaly for M > 1/2, which can be 
observed in the last pane of Fig. |4j The velocity depen- 
dence of the amplitude of oscillations derived from Eqs. 
(15 18) is shown in Fig. [51 It can be seen that, with in- 
creasing velocity, the Friedel oscillations become stronger 
behind the charge (especially for M > 1/2), and weaker 
in front. This is qualitatively illustrated in Figs. [6] and [7] 
Note, however, that no substantial change occurs in the 
Friedel oscillations in the direction perpendicular to the 
motion for any v < v F - 

When the test charge is moving at the Fermi velocity 
(see Figs. [8] and [9]) the stronger Friedel oscillations are 
still present behind the test charge in the far-field, but 
close to the test charge a strong irregular wake field of 
a somewhat different character emerges. In particular, 
this near-field wake is still distinctly present even for very 
small values of r/ (e.g. r\ — 0.02). 



where the first term inside the curly brackets represents 
the semiclassical asymptotic contribution as described in 
Sec. Ill A with H(-y) given by Eq. (Ill; and the sec- 



ond term represents the contribution due to the Kohn 



2. The case of v > v F 

The inclusion of quantum recoil increases the lower 
bound on the phase velocity u/k of undamped longitu- 
dinal oscillations [35] • This can be shown by numerically 
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FIG. 6. The quantity ipr 3 [in units of QAy F /(47reo)], along the direction of motion, in the semiclassical case (blue, non- 
oscillatory), and the case r\ = 0.5 (green, oscillatory), at various velocities of the test charge. 
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FIG. 7. The quantity (pr s [in units of Q\ TF /(4ire )] as a 
function of position for r\ = 0.5, with the test charge at the 
origin and moving at speed v — 0.5«f to the right. The 
black and white lines denote the boundary between positive 
and negative potential for r\ — 0.5 and the semiclassical case 
respectively. 
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FIG. 8. The quantity (pr 2 [in units of QXtf /(4vreo)] as a 
function of position in the case r\ = 0.5, with the test charge 
is located at the origin and moving to the right at the Fermi 
velocity vf- 



solving the equation D(w,k) = for real u and k, where 
D(u>, k) is given by Eq. Q, or its equivalent Eq. Q. The 
minimum phase velocity ^thresh is greater than the Fermi 
velocity, coinciding with it in the semiclassical limit, as 
shown in Fig. [TT] Hence, for test charge velocities in the 
range vp < v < ^thresh) the charge cannot excite un- 
damped oscillations, and the wake depicted in Fig.[T0|is 
formed. This is not as strong as in the semiclassical case 
(it falls off faster than 1/r), and is also of a markedly 
different character, bearing no resemblance to a Mach 



cone. 



For v > Wthresh, on the other hand, all undamped longi- 
tudinal oscillations with phase velocity i>thresh < v^ < v 
are excited. We have not performed calculations for this 
case, but the results of Ref. [48] show a strong oscillatory 
wake, though its asymptotic behavior is not apparent. 
For very large velocities (v 3> vp), the wake formed will 
be identical to the semiclassical wake 1491. 
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FIG. 9. The quantity <pr 2 [in units of QATF/(47reo)] behind 
the test charge for v — Vf- Green is r\ — 0.5; red is r\ = 0.02; 
blue is semiclassical. The test charge is at the right of the 
figure (rn = 0). Inset: an extension of the main graph, for 
r\ = 0.5 (green), and r\ — 0.3 (magenta). 




FIG. 10. The ratio of the potential v?( r ) to the unscreened 
Coulomb potential Q/(4neor), in the case r\ = 0.5, with the 
test charge at the origin and moving at speed v = IAvf to 
the right. 



IV. DISCUSSION 

Let us discuss the applicability of the model. Firstly, 
our model is based on the mean-field approximation 
which is justified at ?] < 1 pS [29l [34] ■ In the limit i] -> 0, 
however, the quantum recoil disappears [34] . Thus the 
fact that our model includes the quantum recoil but does 
not include particle correlations makes it inconsistent in 
a certain sense. Nevertheless, it is widely used in the 
literature, i.e. at small r\ the quantum recoil is assumed 
to be more important than particle correlations. Their 
effects are discussed, e.g., in Refs. [29 | [50U53] . 

Secondly, our model does not include relativistic ef- 
fects. This imposes a lower bound on the parameter r\ 



FIG. 11. The dependence of the threshold velocity Vthresh 
on the parameter r\. In the region marked "undamped" , 
undamped oscillations are excited. In the region marked 
"damped", they are not excited. 



since the latter can be represented as 



n 



a c 
3n vp ' 



(19) 



where c is the speed of light and a = e 2 /(4ireohc) w 1/137 
is the fine structure constant (here we assumed that the 
charge of the particles is equal in its absolute value to 
the elementary charge). The condition vp <C c requires 
r] ^$> yJa/{2>Tr) w 0.03. Hence the limit of rj — > 0, as 
considered in Sec. |IH A[ cannot be taken strictly. Nev- 
ertheless, the results of Sec. |HI A 



are relevant in a cer- 
the limit rj — > corre- 



tain parameter regime. Indeed 
sponds physically to the neglect of the quantum recoil, 
while the quantum recoil provides a contribution given 
approximately by Eq. (13) (for M <C 1 and r 3> Xtf)- 
This is negligible compared to the l/r 3 -term in Eq. (12) 
for 5?7 4 <§; M. (To be exact, the latter inequality is for 
| cos 7 1 ~ 1, while for cos 7 — it changes to 2rj 2 <M.) 



Thus the applicability range of Eq. ( 12 ) due to the above 
restriction r\ > 0.03 is 4 x 10 -6 < 5?T < M < 1 
it is for I cos 7 1 ~ 1, while for 
2 x 10~ 3 < 2ri 2 < M 



(again, 
cos 7 = it changes to 

x iu " « 2rf <C M <C 1). As the velocity increases, 
the effect of the quantum recoil becomes less important 
in front of the charge (such that for v/vp not too close 
to zero the semiclassical and general forms of the poten- 
tial are very similar even for 77 = 0.5, as can be seen in 
Fig.fHJl, and more important behind the charge, especially 
at vjvp > 0.5. 

Thirdly, our model deals with the linearized potential. 
The linearization is justified when the potential energy of 
the interaction of a plasma particle with the test charge 
at the characteristic screening length (i.e. the length at 
which the Coulomb potential and the actual potential 
start to significantly deviate from each other) is much 
smaller than the characteristic kinetic energy of the par- 
ticles in the frame of the test charge. At small veloc- 
ities (v < vp) this condition reads |Qe|/(47reoATF) *C 



Ep which is equivalent to 12tt{\Q / e^if <C 1, while at 
large velocities (v > Vp) the characteristic screening 
length becomes v /ui p and the above condition changes 
to \Qe\ujp/ '(Airegv) -C mv 2 which can be rewritten as 
\2Tr(\Q/e\)rp(vF /v) 3 <C 1. Thus larger velocities imply 
better applicability of the linear theory. 

Finally, there is a number of other effects which are 
not considered in our model but may be important under 
certain conditions. For instance, in our model the shield- 
ing is due to one plasma component (e.g. electrons) only. 
The screening in the presence of both electron and ion 
response is considered in Refs. [SUES] using the quantum 
hydrodynamic model [34] . Furthermore, we did not in- 
clude plasma production and loss processes. In classical 
plasmas they can introduce an additional non-screened 
term and/or a weakly screened term in the potential and 
thus change its long-range behavior, even in the absence 
of the charge motion [56-59 . In quantum plasmas these 
processes (e.g. pair creation and annihilation) can have 
an important effect on the wave dispersion [5D] and thus 
can affect the shielding of a moving charge. 

The potential of a moving test charge under the Lind- 
hard dielectric function has previously been investigated 
in Ref. [JS]- That paper focused on the potential along 
the line of motion and surrounding the first potential 
minimum behind the test charge, for values of r/ approx- 
imately in the range 0.3 - 0.6, while the present pa- 
per is focused on the three-dimensional distribution and 
asymptotic behavior. Since the authors of Ref. [48] con- 
sidered only the "wake potential" , i.e. with the Coulomb 
potential subtracted, the actual behavior of the potential 
is not apparent from their work except when the wake is 
strong enough that the Coulomb potential is small by 
comparison, which occurs only when undamped oscilla- 
tions are excited, i.e. v > ^thresh- We have not studied 
this case in the present paper. 

Our results can be applied, for instance, to investiga- 
tion of bound electron states in the wake fields of ions 
in solids. This problem was studied in Ref. [15] using 
a parabolic fit to the potential around the first poten- 
tial minimum, while a detailed investigation requires an 
accurate knowledge of the whole potential distribution. 
Furthermore, our model can be extended to calculate the 
stopping power or the drag force, similar to classical plas- 
mas HlEiHel. 



V. CONCLUSION 

We have studied the shielding of a moving classical 
test charge in a fully degenerate fermion plasma in both 
the semiclassical case (when the degeneracy of the un- 
perturbed state is the only important quantum effect), 
and the case when quantum recoil is included. 

In the semiclassical case for v < vf , the potential goes 
asymptotically as 1/r 3 , and is repulsive in front of the 
test charge, and attractive behind the test charge and 
perpendicular to the motion; for v close to vp the at- 



traction is especially pronounced perpendicular to the 
motion. For v > vp the test charge excites plasma oscil- 
lations and a strong oscillatory wake is formed inside the 
Mach cone. 

We have also found that the inclusion of quantum recoil 
leads to new effects entirely absent from the semiclassi- 
cal case and the case of a Maxwellian distribution. The 
Friedel oscillations, already present in the screening of 
a static charge, are increased in strength behind a mov- 
ing charge, and for v > vp/2 the asymptotic behavior 
of their amplitude behind the charge changes from I /r 3 
to I/r 25 . Furthermore, the inclusion of quantum recoil 
makes the threshold velocity for excitation of an oscilla- 
tory wake of undamped plasma oscillations larger than 
the Fermi velocity. 

These findings extend the previous results on the 
shielding of a moving charge in a classical plasma to the 
quantum case and can be applied, for instance, to inves- 
tigation of bound electron states in the wake fields of ions 
in solids. 
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Appendix A: Non-dimensional form of the equations 

In order to demonstrate that our choice of free param- 
eters M and r\ determines the potential up to scaling, we 
express the problem in non-dimensional form in terms of 
these parameters. Defining the normalized wavenumber 
and frequency by q = (vf/uj p )Il, w = u/lu p , the dielectric 
function of Eq. Q can be written as 



D(uj,k) = D(w,q) = 1 
x F (w + iv - 



m 2 ,q) 



q 2 2r\q 2 
F (w + i 



vq 2 ,q) 



(Ala) 



where v = v /uj p is an infinitesmal positive number, and 



F(W,q) 



W 



W 2 , [W 
In 



Aq 



\W-q 



(Alb) 



If we now define the normalized variables x = t/Xtf, 
<!> = (47T£oAtf /Q) l P, the non-dimensionalized form of 
Eq. (HI becomes: 



$(x) = 



1 



exp ( zq • x 



2n 2 VSJ q 2 D{Mq-v,q) 



dq, (A2) 



where v is the unit vector pointing in the direction of 
motion of the test charge. 



Appendix B: Derivation of Eqs. (15-181 



larity. For r and v parallel, this is: 



In order to derive the asymptotic contribution to the 
potential from the Kohn anomaly, we first integrate 
Eq. (HI) by parts twice. The surface terms at k = and 
at infinity vanish; in order to avoid having to evaluate 
surface terms at the locations of the Kohn anomaly, we 
take the limit v — >• + only after integrating by parts. 
This leaves: 



ip(r) = - 



1 



47ren 2ir 2 r 2 



/[( 



d_ 

dk 



fc 2 L>(k-v,k) 



exp(ik • r)dk. 

(Bl) 



Assuming now that r is parallel or antiparallel to the 
motion, we can write 



dk^2n dfi 



(B2) 



where /i = k • v. Replacing the integral over k by the 
asymptotic contribution thereto from the locations of the 
Kohn anomaly, k — 2kp(l ± Mfi), we obtain an interme- 
diate expression, which we denote <y9kohn, for the asymp- 
totic contribution to the potential from the Kohn singu- 



Vkohn(r) = i 



Q 



sgn(/i) 



47re 2r 2 A| F J_ x | k\[D{k + v^k + )f 



^d -,=> 



exp(ifc + r^) >dfj,, (B3) 



where k + = 2kp(l + M\fi\). Only the outer surface of 
the Kohn anomaly, i.e. |k| = k + , contributes to this ex- 
pression; the contribution of the inner surface vanishes 
because it involves the integral J^° e lx /(x + iO)dx, which 
evaluates to zero by closing the contour in the upper half- 
plane. For r and v antiparallel, on the other hand, we 
get: 



<Pkohn(r) 



Q 



i 



sgn(/i) 



Aire 2r 2 \ 2 TF J A k 3 _[D(k_vfi,k_)} 



2k F M 



(1 - M 2 ) 



exp(— ik-Vfi) >d/j,, (B4) 



where k- = 2kp(l — M\fi\). This time only the inner 
surface, i.e. |k| = fc_, contributes. 

The asymptotic formulas given in the body of the pa- 
per, Eqs. (pl-pl), follow from Eqs. pi and (|B4| by 



changing the variable of integration to £ = /ifc + (/i) [in 
Eq. (B3)], or £ = /nfc_(/z) [in Eq. (B4)], and considering 
the asymptotic contribution as r — > oo from /i = ±1 
(the endpoints of the integration), and from fi = 
(where the integrand is discontinuous). In addition, the 
stronger ~ i/ r 2 .5 p r j ec Jel oscillations arise from points 
where d£/dfi = 0; this does not occur for r and v par- 
allel, whereas for r and v antiparallel and M > 1/2, it 
occurs at fi = ±1/(2M). 
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